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Abstract 

The states of the three-dimensional quantum harmonic oscillator classify optical 
aberrations of axis-symmetric systems due to the isomorphism between the two mat- 
hematical structures. Cartesian quanta and angular momentum classifications have 
their corresponding aberration classifications. The operation of concatenation of op- 
tical elements introduces a new operation between harmonic oscillator states. 


1 Introduction: Optical Phase Space 

Geometric optics uses the following ‘screen’ coordinates for light rays [1]: 
q = (q X i<ly)’ Coordinates of position on the 2-dimensional screen. The intersection of the ray 
with the screen ranges over 5R 2 . 


P — ( PxiPy ) : Coordinates of momentum (with respect to the same screen). The projecton of 

the ray 3-vector along the ray [of length n, the refractive index of the medium at 
the point] on the plane of the screen is the momentm 2-vector; the component 
normal to the screen is the evolution Hamiltonian (below). See Figure 1. 

Geometric optics has Hamiltonian evolution equations between the canonically conju- 
gate variables p and q. The optical Hamiltonian is 

h = -p x = -yjn(q,z) 2 - |p| 2 . 

(1) 

Plain geometry [2] provides the first Hamilton equation 


dq dh p , , 

iz - a P ~ ° >q ’ 

(2) 

while Snell’s law yields the dynamics of the second Hamilton equation 


dp dh n dn f . 

dz~ aq" p*aq -{p ‘’ 0>p - 

(3) 

We use the Poisson bracket [3] 


,, , v- (*f d 9 df 9g\ 

U ' 9> > y l d qi d Pi d Pi dqJ ’ 

(4) 
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FIG. 1 . The coordinates (p, q) of optical phase space in geometric optics. The 
z-axis is traditionally the optical axis. 

for functions /, g of p and q, and the corresponding Lie-Poisson operators {/, o} as generators 
of Lie algebras [4]. This space is subject to Hamiltonian evolution equations and constitutes 
optical phase space, subject to the same local symplectic structure as the well known mechanical 
phase space. Optical phase space differs from the mechanical one globally , however, in that its 
momentum ranges over a compact domain: (two) disks (forward and backward rays) of radius 
n(q), the refractive index. 

2 Examples 

Free flight in homogeneous medium (n constant) yields the solutions 

q(z) = exp z{pz , o} q = q + £ . P , p(z) = exp z{p s , o} p = p. (5) 

V” - IpI 

This is shown in Figure 2. 

One example we may use to distinguish mechanical from optical phase space trans- 
formations pertains the mechanical oscillator versus — or vis-d-vis [5] — li ght in a fib er whose 

refractive index is a function of the radius |q| to the optical axis: n(q) = — HqP ■ We ca ^ 
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FIG. 2. The deformation of phase space due to free flight between two screens. 
Since the optical transformation is canonical, the area elements are preserved, and 
so is the allowed strip of phase space, |p| < n. [Figure by G. KROTZSCH, IIMAS- 
UNAM.] 


this the elliptic index-profile fiber. The evolution Hamiltonian is h = — •y/nf - (IpP + Hql 2 )- 
For every z phase space is the patch (|p| 2 + i/|q| 2 ) < n 2 . Figure 3 shows the resulting evolution 
of phase space, compared with that of the mechanical harmonic oscillator. 


In a homogeneous optical medium, phase space is a strip |p| < n. The Fourier 
transform is a fundamental ingredient for much of coherent-state (paraxial) optics; yet, the 
Fourier transformation (p *-*■ — q, q >— ► p) is not an invariance transformation of this space. 
Nevertheless, one may define a truly optical Fourier transform, that in one dimension has the 
form 


. _ d-p 2 )^ , 


[i + (i-pWl 8/2 
v/w 


(6) 


that respects the strip of phase space, paraxially rotates by jir (as the usual Fourier transform), 
and is nonlinear symplectic. In Figure 4 we show the optical Fourier transform for a quadrant 
in the phase space strip. 


3 Linear Transformations of Phase Space 

Although mechanics and optics phase spaces differ in their global properties, their paraxial 
correspondence motivates that we Taylor-expand all expressions into series where we may trun- 
cate the series to some aberration order in the powers of the phase space variables. We will 
thus work with polynomials and so we may ignore the range restriction. In fact, we thus substi- 
tute the structure of Euclidean-based Arc optics by the metaxial Heisenberg- Weyl Lie-algebraic 
structure suitable for aberration expansions in powers of the phase space coordinates. 
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FIG. 3. Left : evolution of mechanical phase space under harmonic oscillator 
time evolution, for one space and one momentum dimension. Phase space rotates 
rigidly. Right: evolution of optical phase space under z-translation. The phase 
space disk rotates differentially and resembles the harmonic oscillator only for small 
p and q, i.e., in the paraxial optical regime. The circular (in general, elliptic) region 
is optical phase space: coordinates axe meaningless beyond the edge. [Figure by G . 
KROTZSCH, IIMAS-UNAM.] 


In effect, we may classify the symplectomorphisms of phase space by the Lie-Poisson 
generators of the finite transformations. 

Translations of phase space are generated by linear functions of (p, q) through expo- 
nentiated Lie-Poisson operators: 


exp(a.{p,o» (*) = („!,), exp(b.{q,o}) (P) = ( p + b ), ( 7 ) 


as shown in Figure 5. 

Linear canonical transformations of phase space are generated by the quadratic func- 
tions PiPj, Piqj , <7 t 9y, i = x, y. These functions close into an sp(4, 5ft) algebra under Poisson 
brackets. Further, optical systems that have a common axis of rotational (and inversion) sym- 
metry are generated by linear combinations of p 2 = |p| 2 , p*q, and g 2 = |q| 2 , that close into an 
sp(2, 5ft) = so(2, 1) = sl(2, 5ft) Lie algebra. The corresponding group transformations of phase 
space are well known and shown in Figure 8. The quadratic polynomial p 2 + q* generates rigid 
rotations, as is well known. 
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FIG. 4. Left : a quadrant in the phase space strip p < 1. The usual (mecha- 

nical) Fourier transform rotates this patch by j* in the counterclockwise direction. 
Right: the optical Fourier transform maps points in the strip on points in the same; 
paraxially it rotates the plane by jt, and is globally symplectic: it conserves area 
elements. [Figure by G . KROTZSCH, HMAS-UNAM.] 

4 General Nonlinear Transformations 

Nonlinear symplectic transformations of phase space are generated by polynomials in the com- 
ponents of (p, q) of degree higher than second. Again, the transformations of axis-symmetric 
optical systems leads us to concentrate on polynomials of three variables 

£+=p 2 /\/2, £o=P-q. £- = 9 2 /V2. (8a) 

— * 

These variables may be placed into a vector £ of cartesian components 

£i = i(p 2 -9 2 )> £2 = j(p 2 + « 2 )’ 6 = p*q- (8&) 

They close into an sp(2,9t) algebra, as we noted before. The Casimir of this algebra is minus 
the squared radius of a sphere, 

(pxq) 2 = pV - (p . q) 2 = -£j - £ 2 - £ 2 = -|£ | 2 , (9) 

and is the well-known Petzval invariant of geometric optics. 

In representing the action of optical elements by means of Lie-Poisson transformations 
generated by polynomials, we are aided by the following theorem: 
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FIG. 5. Translations of phase space are generated by linear functions of the 
phase space observables. Left : exp{p, 0 } translates in q. Right: exp{p, 0 } translates 
in p. The latter does not leave the optical phase space invariant, but serves as a 
first-order approximation to translation and rotation of the screen. [Figure by G . 
KROTZSCH, IIMAS-UNAM.) 




FIG. 6 . Linear transformations of phase space are generated by quadratic poly- 
nomials in the phase space variables. In one position and one momentum coordinate, 
the transformations shown above correspond to exp{p 2 ,o}, exp{pg, o}, exp{<z 2 ,o}, 
and respectively. [Figure by G. KROTZSCH, IIMAS-UNAM.] 


Theorem (Dragt & Finn) [6]: Canonical transformations .M leaving the origin invariant, 
(i.e., excluding translations) can be approximated by a truncated product series of Lie trans- 
formations 

>( = ••• exp{/ 5 ,o} exp{/ 4 ,o) exp{/ s ,o} exp{/ 2 ,o}, (10) 

where /jfc(p,q) are homogeneous polynomials of degree k in the components of p and q. 
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Correspondingly, for axis-symmetric systems, the approximation is written as 

>! = ••• exp{a 4 ,o} exp{a s ,o} exp{a 2 ,o} expfa^o}, (11) 

with a*(£) polynomials of degree k in the components of £ as defined in Eqs. (8); k is called 
the rank of the phase space generating monomials. 

The rightmost factor in (10) and (11) is the subgroup of linear transformations, and 
factors to the left of it are nonlinear, called generically aberrations in optics. Since the degrees 
of polynomials in Poisson brackets satisfy {/y,/*} = /y+*_ 2 , the ranks of the sp(2, &)-based 
polynomials satisfy 

i a j> a k} = °y+Jfc-i* (12) 

it follows that for each rank k, the polynomial set a* is a (reducible) ideal under linear (paraxial) 
transformations. This is the aberration ideal of order A — 2k — 1. Axis-symmetric optical 
systems are described by a paraxial approximation modified by aberrations of orders 3, 5, 7, 

. . . , generated by a 2 , aj, a 4 , . . . , respectively. 

There are legitimate questions about the convergence of the product series and global 
properties of the group of nonlinear symplectic transformations. In any case, by taking functions 
of phase space modulo polynomials of degree higher than the aberration order, we may construct 
a well-defined finite-parameter group of transformations of phase space truncated to that order. 
Its best parametrization and, especially, the product law must be discovered — once and for 
all. Thus we construct the A-th order aberration groups. The theory developed in optical 
aberrations serves as well in higher approximation theory. 


5 The Monomial Classification of Aberrations 


It should be quite evident by now that we may classify aberrations of axis-symmetric systems 
through proposing complete bases of polynomial functions a*(£) of degree 2k in phase space. 
(This may also be used to classify non-axis-symmetric aberrations as broken symmetry [7]. We 
have a sp(2, dt) Lie algebraic graded covering structure with Poisson (Berezin) brackets between 
polynomial functions of three variables. 

The monomial basis is 

^k+,ko,k- = const x £_^. + £ q ° £_ , kj = 0, 1, ... . (13) 

These monomials have rank k = k+ + A;q + k~, and are classified as the harmonic oscillator 
Cartesian basis states. We may examine the action generated by exponentials of these mono- 
mials up to the first Taylor term, 


«P{A/ i+ ,t 0 , t .,°}(^) = ( 


P + *oAf i+iio _ Iit _p + 2*-Af t+ito4 .. 1 q + --\ 
q - 2A:+Afjfc + _ 1( j fcoi j k _p - M4 +) Jfc 0 -l ) Jfc_q+ • • • / 


(14) 


In this Dragt [8] recognized the traditional third-order Seidel aberrations, generated by and 
called 


(p 2 ) 2 SPHERICAL ABERRATION, 
p 2 p • q CIRCULAR COMA, 
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FIG. 7. The third order aberration sextuplet classified in £+, £o> a^d coor- 
dinates, and the traditional names. S: spherical aberration, C: coma, A: astigmatism, 
F: curvature of field, D: distorsion, and P: pocus. 


(p • q ) 2 ASTIGMATISM, 

p 2 q 2 CURVATURE OF FIELD, 
p • q q 2 DISTORSION, 

(g 2 ) 2 POCUS. 

Actually, pocus was excluded from Dragt’s list because it does not produce any change 
in the screen images; rather it only changes the directions of ray arrival; it has a p-unfocusing 
action that we have playfully called pocus [9]. Yet it is an aberration of phase space on par 
with the other five, and the Fourier transform of spherical aberration. In Figure 7 we show the 
familiar harmonic oscillator states with two quanta and its corresponding aberrations of rank 
two in the monomial basis. 

6 The Symplectic Classification of Aberrations 

The harmonic oscillator eigenfunction structure naturally suggests another basis, following the 
number of quanta, total angular momentum, and one projection. 

The symplectic harmonic basis is [1], [9], [10], [11] 

*4,(0 = const x [({) J | ( ‘- y)/! V,>( f ), (15) 

where (£ ) 2 is the Petzval and }/j im (£) the solid spherical harmonic of angular momentum j 

and projection m. The total degree of *£•«(£) in is k —the rank. The index j we may call 
the symplectic spin , and m the Seidel index [12]. The third-order aberration sextuplet thus 
reduces, under linear transformations, into a spin-2 quintuplet and a spin-0 singlet (that is in 
effect the Petzval itself). Only the m — 0 member of the quintuplet and the singlet mix the 
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FIG. 8. The aberration multiplets to seventh order classified into symplectic 
harmonics by rank, spin and Seidel index m. The symbols are 0: sphericalsabe- 
rration, O: circular coma, O, >: elliptic comas, CD, txi: curvatisms, |, cx: 
astigmatures, •*- k : distorsion, and •: pocus [12]. 


curvature of field and astigmatism monomials into aberrations that we could call curvatism (in 
the quintuplet) and astigmature (the singlet, invariant under paraxial transformations). The 
scheme appears in Figure 8; this is but the harmonic oscillator t-m spectrum from the k = 2 
level (s-d shell) up, seen sidewise. 

The constant in front of the definition of *Xm in Eq. (15) has been chosen to avoid square 
roots, because they are time-consuming for symbolic manipulation. Thus we have defined them 
starting from the normalization of the highest weight = (p 2 )*, and their lowering through 

{ Ifl 2 . °} = (rn + j) They we 




An(2j -)- 1) (j + m)\ (j - m)l 

N (2 j - 1)H 

Qjjj ~ rn) 

2 m / z (2j - 1)!! „ 


yL( o 


2 m / 2 f2» - I’ll! ,2 n (m + n)! (j-m- 2n)! n! 


We give the symplectic harmonics only for m > 0,, because 

o,e-) = *x£(e-,e 0 ,e+). 


(16) 


(17) 


The first two ranks ( k = 0, 1) correspond to the scalar and the sp(2, 5i) generator 
functions: 

%° = 1; l Xh = V2£ ±t 'X} = fa (18a) 
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The Jfc = 2 basis functions are the generators of third-order aberrations: 

2 4 = nl, 

2 X} = \/2£+ fo, 

% 2 = i(f+e- + f 0 2 ). 2 r 0 ° = 2 «+«. - «g. 


For Jfc = 3 we have the generators of fifth-order aberrations: 

s rf = 2v^d. 
s 4 = nl(o, 

5 4 = |v5(el€- + «+«?), s 4 = J 4° *4. 
*X$ = i(3€+6)C- + «J). = 2 4° *4- 


For Jfc = 4 we have those of seventh-order aberrations: 

‘xi = 4ft, 

*X , 4 = 272ft f 0 , 

‘4 - f (fif- + 3f 2 ft), 4 4 = 2 4 2 4, 

44 = ?72(3ftfof- + 2f+fo). 4 r? = 2 4 a 4. 

4 4 = &(3ft ft + i2f+f 2 f- + 2ft), 4 4 = 2 4 2 4, %° 


(W 


(186) 


(18c) 


(18«0 


7 Spot Diagrams of Harmonic Oscillator States 

The epot diagram of a tramformation M : (p, q) >-» (p'(p, qW(p,q)) i» > projection of S 4 
phase space on the screen plane, that pictures q'(p,q), the image of a pencil of rays (range 
of p) diverging from a fixed object point q. If we let p mark & polar coordinate grid around 
the optical axis, we obtain the spots of Figures 9.1-9.5 [12]. These axe the new “faces” of the 
harmonic oscillator states that we present in aberration optics. 


8 Characterization of Optical Elements 

Optical elements may be characterized by 

£{A 4 ,Aj, A 2 ,M> = exp{A 4 ,o} exp{i4 s ,o} exp{A 2 ,o} exp{A lt o}, (19a) 

with the coefficients their linear action M = d)’ = coefficients Akj tm 

of its aberration polynomials 

Ajfc= E° E Akj,m k xLii). m 

j=k,k- 2 ,... m=j, 
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FIG. 9.1. Spherical aberrations © generated by (at left we indicate the 
kjm classification of the spot diagram). 


Some of the most common optical elements thus represented are the following: 

Free propagation: The paraxial part is F(a) = ^ \j n The (spherical aberration) 

coefficients are: 

fw = -*/(8* S ), /s. 3,3 = -V(16n 5 ), / 4 Z , 4,4 = -5z/(128n 7 ). 


Elliptic-profile fiber free propagation in medium n = 


is treated in [13]. The 
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FIG. 9.2. Circular comas O generated by k X k _ v 
[ BE • HE \ 

paraxial part is F n(g) (z) = [ ^ aberTati ° n P ol y nomials are: 

F 2 = -z(|p| 2 + M 2 |q| 2 ) 2 /8^, 

Fi = - 2 (!pI 2 + M 2 |q| 2 ) s / 16l/5 » 

F 4 = -5z(|p| 2 + ^ 2 |q| 2 ) 4 /128i/ 7 . 

The root transformation [14] indicated by (that is the root of refraction see below) 

in medium n associated to the surface S(g 2 ) = &9 2 + + **"• ^e P 81831 * 8 ! P 81 ^ 

r c _ ( 1 _2n?2 ^ . The aberration coefficients R k .• m , arranged as row multi-vectors with 

n >* \0 1 / 

components numbered by descending values of j and m , are 

2 r = {{0,0, -ft/(2n), 0, n? 4 >, {-? 2 /(3n)}}, 
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FIG. 9.3. The curvatism-astigmature m-degenerate pairs CD — |, generated 
by k Xjf_ 2 and k X k Z%- 


3 r = {{0,0, -f 2 /(8n s ), -<T^/(2n 2 ), -f 4 /(2n), 2&f 4 , nfc}, {-? 2 /(10n 3 ), -?£/( 5n 2 ), -2/5 ? 4 /n}}, 

4 r = {{O,O,-fc/(10n 5 ),-^/(4n 4 ),-f 4 /(8n 3 ) - 5/0?|/n 3 , -2f 2 ? 4 /» 2 , -fc/(2n) + 
16/3ffc 4 /n,ef 2 fc,8/3nf 2 £? + nft}, {~3/50f 2 /n 5 , -^/(7n 4 ), -f 4 /(7n 3 ) - 
2/7?!/« S , -8/7<^ 4 /n 2 , -3/7fc/n - 16/21 ?2 2 <r 4 /n}, {-<r 4 /(15n 3 )}}. 

The refracting surface 5 between medium n and medium m is given by the factorization 
theorem [14], [15] 

Sfi.mjS = %n,S %m,S’ (20 


The linear part is S = 2m?2 1 2nf2 ^ an d the aberrations 
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FIG. 9.4. The m-degenerate pairs x — X* 


2 s = {{0,0, ft /(2m) - ft/(2n),2 nft 2 /m - 2ft 2 , -mft + 2 + nft ~ 4nft 3 + 2 n 2 ft 3 /m>, {ft/(3m) - 
ft/(3n)}}, 

3 s = {{0,0,ft/(8m 3 ) - ft/(8« 3 ),nft 2 /m 3 - ft 2 /(2m 2 ) - ft 2 /(2n 2 ),3n 2 ft 3 /m 3 -3nft s /m 2 + fc/(2m) + 
2f 3 /m — ^4/ (2rt) — 2f| /n,— 2mftft/n + 4mf2/ n — 2ftft + 4 — 0n ft /m + 4nftft/m + 
4nft 4 /m - 6ft 4 , -mft + 6mft 2 ft - 2mft 5 + nfc - 12nf 2 fc + 4nft 5 + 2n 4 f|/m s -4n 5 ft 5 /m 2 + 
6n 2 ft 2 ft/m},{ft/(10m s ) - ft/(10n 3 ),2/5nft 2 /m 3 - f|/(5m 2 ) - fj/(5n 2 ),2/5n 2 ft 3 /m 3 - 
2/5nf 3 /m 2 + 2/5 ft/m - 2/5?J/m - 2/5ft/n + 2/5ft/n», 

4 s = {{0,0, ft /(16m 5 ) - ft/(16n 5 ),3/4nft 2 /m 5 - ft 2 /(4m 4 ) - ft 2 /(4m 3 n) - ft 2 /(4n 4 ),15/4n 2 ft 3 /m 5 - 

5/2nft 3 /m 4 +ft/(8m 3 )-5/12ft 3 /m s -ft/(8n s )-5/6ft 3 /n s ,10n 3 ^/m 5 -10n 2 f|/m 4 +3nftft/m + 
8/3nc2/m 3 -2ftft/m 2 -4f£/m 2 +ftft/(mn)+4^/(mn)-2ftft/n 2 -8/3^/n 2 ,-10/3mftft/n + 
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FIG. 9.5. Elliptical comas and their m-degenerate partners ( > — >. 


32/3mf|/n 2 + 15n 4 fj/m 5 — 20n 8 f|/m 4 + 15r» 2 f|f4/m s + en 2 ^/"* 5 — 62/3nf 2 f4/"i 2 - 
20/3 nfj/m 2 + fc/(2 m) + 19f|f 4 /m + llf|/m - ?<j/(2n) - 8^f 4 /n - 16s$/n,-6m$2&/n + 
116/3m£|f 4 /n - 52/3 mfj/n + 4f 2 fc + 12 n 5 fj/m 5 - 20 n 4 fj/m 4 + 28n s f|? 4 /m 5 - 4/3 n s f|/m s - 
176/3»t 2 f|f 4 /m 2 + 64/3n 2 f®/»Ti 2 + 2nf2fc/m + 72n$2(i/m — 36 rt^/m + 4n{ ,2 /m — 80f|f4 + 
124/3f| - 4^,8mf 2 ? 4 - mfg + + 86/3mf£f 4 - 124/3 mfj - 24n?2$ 4 + nft - 4n^ fl - 

40n^f 4 + 200/ 3rtfJ + 4 »®$J /m 5 - 8n 5 /m 4 + 18n 4 ^ 4 /m 2 - 20/3» 4 fJ /m s - 152/3n S f 2 4 ?4/m 2 + 
112/3n 2 $« /m 2 + 40/3n 2 tf2f 2 /m + 2n 2 f 2 Gb/th + 164/3n 2 f<2 £4/ fn — 188/3n 2 fJ /w» + 8/3m 2 £ 2 £ 4 /n — 
32/3m 2 f2?4/» + 32/3m 2 fJ /^}»{3/56f2/m^ — 3/56^/n®,3/7nf 2 /m® — £ 2 /(7m 4 ) — f 2 /(^ m ^ n ) — 
f2/(7n 4 ),9/7n 2 f2 5 /"* 5 -0/7n^/m 4 +f 4 /(7m 3 )-fJ/(7m s )-f 4 /(7» s )-2/7f2V« 3 42/7n s ^/m 5 - 
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12/7n 2 ?|/m 4 + 12/7 n? 2 ? 4 /m 3 - 8/7 n?|/m 3 - 8/7? 2 ? 4 /™ 2 + 12/7?£/m 2 + 4/7? 2 ? 4 /(mn) ~ 
12/7f|/(mn) - 8/7? 2 ? 4 /n 2 + 8/7?£/n 2 ,16/21m?J? 4 /n 2 - 32/21 m?|/n 2 + 0/7n 4 f|/m 5 - 
8/7 n 3 ?f/m 4 + 20/7n 2 ?f ? 4 /m 3 - 20/7 n 2 ?|/m 3 - 64/21 n?|? 4 /m 2 + 104/21 n?f /m 2 + 3/7 ft/m - 
12/7? 2 ? 4 /m - 18/7? |/m - 3/7? e /n + 8/7? 2 2 ? 4 /n + 16/7?|/n},{? 4 /(15m s ) - ? 4 /(15n 3 )}}. 

The curved mirror transformation by 5 can be found from At n> 5 = £n,S^-n,S [^®1* The 
refracting surface between two elliptic-profile fibers has been calculated in Ref. [13]. 

9 Multiplication Law in the Aberration Group 

For aberration orders 1, 3, 5, and 7 (ranks fc = 1, 2, 3, and 4), the dimensionality of the 
basis is, respectively: the 3 generators of sp(2, 2R), the 0 third-order aberrations (separated 
into a quintuplet and a singlet), the 10 fifth-order aberrations (a septuplet and a triplet), 
and the 15 seventh-order aberrations (divided into a nonuplet, quintuplet and singlet). The 
number of parameters of the corresponding aberration group elements in Eqs. (11) and (19) 
thus accumulates to: 

3 SP(2,D?) (ABERRATION ORDER 1) 

9 LINEAR + ABERRATION ORDER 3 
19 UP TO ABERRATION ORDER 5 
34 UP TO ABERRATION ORDER 7 

If we indicate the seventh-order pure aberration group elements by the coefficients of 
the polynomials in 

5{A 4 ,Aj, A 2 ,l} = exp{Ai,o} exp{Aj,o} exp{A 2 ,o}, (21) 

the central problem is to find the multiplication table involving the 31 up-to-seventh order 
aberration coefficients in the product 

£{C 4 ,Cj,C 2 , 1} = £{A 4 , A}, A 2 ,1}5{B 4 ,Bj,B 2 ,1}, i.e. C = A # B. (22) 

To find explicitly the goto operation # between the individual coefficients, we may use 
Baker-Campbell-Hausdorff relations. Order 3 is abelian; order 5 is the practical limit for 
hand calculations, and order 7 is definitely nontrivial and needs symbolic computation [12], 
[16]. The pure aberration group composition law was calculated once and for all, to find the 
composition of aberrations. It is: 

Aberration order 3: 

^2 ,j,m ~ J = 2, 0, m — J,J — 1, • . . , ~j< 

Aberration order 5: 

<?3,3,3 = 2A 2)2( iB 2|2)2 - 2A 2(2i2 B 2)2 ,i + -^3,3,3 + -Bs.S.Si 
<?3,3,2 = 4A 2 j2i 0 B 2i2)2 - 4A 2 i 2|2 B 2i2i0 + ^3,3,2 + #3,3, 2» 
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£s,S,l = 642,2,-102,2,2 + 2^,2,002,2,1 ~ 2^2,2,102,2,0 “ 6^2,2,202,2,-1 + -^2,3,1 + 03,3,1, 

C'S.S.O = 8^2,2,-202,2,2 + 442,2,-102,2,1 ~ 442,2,102,2,-1 ~ 84 2| 2,202,2,-2 + -^3,3,0 + 03,3,0, 
Cs.S.-l = 6^2,2,-202,2,1 + 2^2,2,-102,2,0 ~ 2^2,2,002,2,-1 ~ 64 2 ,2,102,2,-2 + 4j, 3,-1 + 03, 3,-1, 
£3,3-2 = 4^2,2,-202,2,0 “ 4^2,2,002,2,-2 + ^3,3, -2 + 03, 3, -2, 

£s,3,-3 = 2^2, 2, -202, 2,-1 - 242,2,-102,2,-2 + ^3,3, -3 + 03,3,-3; 

C's.l.l = 4/54 2 ,2,-102,2,2 ~ 2/542,2,002,2,1 + 2/542,2,102,2,0 “ 4/542,2,202,2,-1 + ^3,1,1 + 03,1,1. 

C'S.I.O = 16/542,2,-202,2,2 ~ 2/542,2,-102,2,1 + 2/542,2,102,2,-1 “ 16/542,2,202,2,-2 + 4s,i, 0 + 
03,1,0. 

£s,l,-l = 4/542,2,-202,2,1 “ 2/542,2,-102,2,0 + 2/542,2,002,2,-1 ~ 4/542,2,102,2,-2 + -^3, 1,-1 + 

*».W 


Aberration order 7: 

Q.4,4 = 842,2,102,2,2 + 32/34? 2 , 202,2,0 " 16/342,2,002,2,2 ~ 32/342, 2, o42, 2, 202,2, 2 “ 

84 2 .2.142.2.202.2.1 +44 2 ,2,102,2,102,2,2 +64 2 ,2,103,S,3 -442,2,202,2,1 + 16/342,2,202,2,002,2,2 ~ 

442,2,203,3,2 + 44,4,4 + 04,4,4. 

^4,4,3 = 324^2,202,2,-1 “ 1642,2,-102,2,2 ~ 3242, 2, -l42,2, 202, 2, 2 + 64/342,2,o42,2,102,2,2 “ 

80/342,2,0-42,2,202,2,1 - 8/342,2,002,2,102,2,2 + 1242,2,003,3,3 + 16/342, 2, M 2 , 2, 202, 2,0 + 
40/342,2,102,2,002,2,2 + 242,2,103,3,2 + 1642,2,202,2,-102,2,2 ~ 32/342,2,202,2,002,2,1 ~ 

842.2.203.3.1 + 4 4i4> s + 04,4,3. 

£4,4,2 = 64/34^2 n02, 2, 2 + 8/34? 21 i?2, 2,0 + 644? 2,202,2,-2 “ 3242,2,-202,2,2 - 

6442.2, -242,2,2-02,2,2 + 1642,2, -l42, 2, 102, 2, 2 ~ 5o42,2,-l42,2,202,2,l ~ 2042,2,-102,2,102,2,2 + 

1842.2, -103,3,3 ~ 4/342,2,002,2,1 “ 8/342, 2, o42, 2, 102, 2,1 “ 64/342, 2, o42, 2, 202, 2,0 + 
32/342,2,002,2,002,2,2 + 842,2,003,3,2 + 4042, 2 , l42, 2, 202, 2,-1 + 2842,2,102,2,-102,2,2 + 
4/342,2,102,2,002,2,1 ~ 242,2,103,3,1 “ 32/342,2,202,2,0 + 3242,2,202,2,-202,2,2 ~ 

842.2.202.2, -102,2,1 ~ 1242,2,203,3,0 + -^4,4,2 + 04,4,2. 

^4,4,1 = 16/34? 2,002,2,1 + 164? 2,102,2,-1 ~ 9642,2,-242,2,202,2,1 ~ 4842,2,-202,2,102,2,2 + 

2442.2, -203,3,3 ~ 842,2,-l0?,2,l + 160/342, 2, -l42, 2, 002, 2, 2 ~ 1642,2,-l42,2,102,2,l “ 
176/342,2,-1 42,2,2 02,2,0 “ 8/342,2,-102,2,002,2,2 + 1442,2,-103,3,2 ~ 16/342,2,042,2,102,2,0 + 
16/342, 2, 042, 2, 202, 2,-1 + 88/342,2,002,2,-102,2,2 + 8/342,2,002,2,002,2,1 + 442,2,003,3,1 ~ 
8/342,2,10?,2,O + 9642, 2, l42, 2, 202, 2, -2 + 4842,2,102,2,-202,2,2 + 842,2,102,2,-102,2,1 ~ 
642,2,103,3,0 ~ 80/342,2,202,2,-102,2,0 ~ 1642,2,203,3,-1 + 44,4,1 + 04,4,1. 

^4,4,0 = 404? ,2,-102,2,2 + 404? 2 , 1 .B 2 , 2, -2 — 2042,2, _202,2,i + 160/342,2,-242,2,002,2,2 ~ 

4042.2, -242,2,102,2,1 - 320/342,2,-242,2,202,2,0 “ 80/342,2,-202,2,002,2,2 + 2042,2,-203,3,2 + 

40/342, 2, -l42, 2, 002, 2,1 - 80/342, 2, -l42, 2, 102, 2,0 “ 4042, 2, -l42, 2, 202, 2,-1 + 

2042.2, -102,2,-102,2,2 ~ 20/342,2,-102,2,002,2,1 + 1042,2,-103,3,1 + 40/342, 2, o42, 2, 102, 2,-1 + 
160/342, 2 , o42, 2, 202, 2, -2 + 160/342,2,002,2,-202,2,2 + 40/342,2,002,2,-102,2,1 + 

2042.2.102.2, -202,2,1 “ 20/342,2,102,2,-102,2,0 ~ 1042,2,103,3,-1 ~ 2042,2,202,2,-1 “ 

80/342,2,202,2,-202,2,0 - 2042,2,203,3,-2 + -^4,4,0 + 04,4,0. 

£4,4,-l = 164? ,2,-102, 2,1 + 16/34? 2(O 02, 2,-1+9642,2,-242,2,-102,2,2 + 16/342,2,-242,2,002,2,1- 
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176/3^2, 2-2-42, 3,0 - 96X2,2,-2-42,2,232,2,-1 " 80 / 3 ^, 2,-232,2,032,2,1 + 16^2,2,-233,3,1 “ 
8/3^2,2,-l3j ( 2 j0 - 16/3^2,2,-1^2,2,032,2,0 “ 16X2,2,-l-42,2,l32,2,-l + 48X2,2,-l3 2 ,2,-232,2,2 + 
8A2,2,-l32,2,-l32,2,l + 6X2,2,-133,3,0 + 160/3^2,2,0^2,2,l32,2,-2 + 88/3^2,2,032,2,-232,2,1 + 
8/3^2,2,032,2,— 13 2 , 2,0 - 4 ^2,2,o3s, 3,-1 “ 8^2,2,l3^2,_x ~ 8/3A 2t 2,l B 2,2, -2 B 2, 2,0 ~ 
14^2,2,133,3,-2 ~ 48^2,2,232,2,-232,2,-1 “ 24A 2> 2,2 B 3,i,-S + ^4, 4,-1 + 3 4 ,4 _i, 

C4,4,-2 = 64A^ 2,-232,2,2 + 8 /3A^ 3 _ 1 5 2 ,2,0 + 84jZA\ 2 0 32,2,-2 - 32/3^2,2,-231,2,0 + 

40X 2 , 2 , -2^2,2,-132,2,1 ~ 64/3^2,2,-2^2,2,032,2,0 “ 56X2,2,-2^2,2,132,2,-1 ~ 
64^2,2,-2^2,2,232,2,-2 + 32^2,2,-232,2,-232,2,2 ~ 8 X 2 , 2 , -232, 2 , -l32, 2,1 + 12X2, 2 , -23s, 3,0 ~ 
8/3A 2 , 2 , -1^2,2,032,2,-1 + 16^2,2,— 1-^2, 2, l3 2 , 2, —2 + 28X2,2,-l3 2 ,2,-232,2,l + 

4/3^42,2,— 132, 2,-1 32, 2,0 + 2X2,2,-l3s,3,-l — 4/3X2, 2, 03^,2, -l + 82/ZA2 t 2,0 B 2,2,-2 B 2,2,0 ~ 
8-^2,2,o33,3,-2 ~ 20^2,2,132, 2, -232, 2,-1 — ^A2 t 2,l B 3,3,-3 — 32X2,2,231,2,-2 + ^4, 4, -2 + 34,4, -2* 

<74,4-3 = 32A^ 2 ,- 2 32,2,l + 16/3^2,2,-2^2,2,-132,2,0 ~ 80/3^2, 2, -2^2, 2, o3 2 , 2,-1 ~ 

32i42,2,-2- / l2,2,l32,2,-2 + 16^2,2,-232,2,-232,2,1 - 32/3X2, 2, -232, 2, -l32, 2,0 + 

8A 2) 2-2 B 3,3,-l + 64/3^2, 2 , -1^2,2,032,2,-2 + 40/ZA 2t 2,-l B 2,2,-i B W ~ 2 A 2,2,-l B 3,3,-2 ~ 
8/3^2,2,032, 2, -232, 2,-1 _ 1 2 A2,2,0 B 3,3,-3 ~ ^^2,2,l B 2,2,-2 + ^4, 4, -3 + B 4, 4,-3. 

C’4,4,-4 = 32/3^,2,-232,2,0 + 8^^2,-l 5 2,2,-2 ~ 4X2,2,-232,2-1 ~ 8A2,2,-2A 2 ,2,-l B 2,2,-l ~ 

32/3i4 2 , 2, -2^2,2,032,2,-2 + 16/ 3^2,2,-232,2,-232,2,0 + 4X2,2,-23s,3,-2 + 

4^42,2,-132,2,-2-32,2,-1 “ 6X2,2,- l3s, 3, -3 - 16/3-42,2,031,2,-2 + -^4, 4, -4 + B 4,4,-4\ 

<74,2,2 = ~64/21-42,2,o32,2,2 ~ 8/21^2,2,132,2,0 + 256/21^2 2, 232, 2, -2 - 128/21.42, 2, - 2 B 2 , 2,2 ~ 
250/21.42,2,-2^2,2,232,2,2 + 170/21X2, 2, -1-42, 2, l32, 2, 2 “ 16/3^2,2,-1^2,2,232,2,1 + 
32/21-42,2,-l32,2,l32,2,2 + 24/7X2,2,- 1 3s, 3, 3 + ^/ 2l ^2,2fi B 2,2,l + 8/21X2,2, 0-42, 2, l32, 2,1 + 
64/21-42,2,0^2,2,232,2,0 ~ 32/21^2,2,o32,2,o32,2,2 “ &/7 A 2t 2,0 B 3,3,2 “ 64/21^2,2,1^2,2,232,2,-1 + 
8/3-42,2,132,2,- 132, 2, 2 “ 4/2lA 2t 2,l B 2,2,0 B 2,2,l + 2A 2 t 2,l B 3,\,l + 24/35^2,2,l33,S,l + 
32/21^2,2,23^2,0 + 128/21^2,2,232,2,-232,2,2 “ 88/21X2, 2, 232, 2, -l32, 2,1 ~ 4X2,2,233,1,0 “ 
24/35X2,2,23s,s,0 + X4,2,2 + 34,2,2 » 

<74,2,1 = -16/7 X^ 2,o32,2,1 - 88/21X2 2, i 32, 2,-1 + 64/3X2, 2, -2-42, 2, l32, 2, 2 - 
704/21X2,2,-2-^2,2,232,2,1 — 128/21X2, 2, -232, 2, l32, 2, 2 + 90/7X2, 2, -2 3j, 3, 3 + 

44/21X2, 2 , -i3 3 , 2,1 - 32/21X2,2,-1X2,2,032,2,2 + 88/21X2,2,-1X2,2,132,3,1 “ 
32/21X2,2,-lX2,2,232,2,0 ~ 32/21X2,2,— 132, 2, 032, 2, 2 + 16/7X2, 2, 0-42, 2, l32, 2,0 + 
64/21X2,2,0-42,2,232,2,-1 + 16/21X2,2,o32,2,-l32,2,2 ~ 8/7 A2 t 2fi B 2,2fi B W + 4X 2 ,2,o3s,l,l - 
32/35X2, 2, o3s, 3,1 + 8/7X2, 2, l3 3 , 2,0 + 256/2lX2,2,lX2,2,232,2,-2 + 352/21X2, 2, l32, 2, -232, 2,2 “ 
44/21X2,2,l32,2,-l32,2,l “ 2X2,2,l33,i,o + 48/35X2,2,l33,3,o - 32/3X2,2,232,2,-232,2,1 + 
16/21X2,2, 232, 2, -l32, 2,0 - 8X2,2,233,1,-1 - 96/35X2,2,233,8,-1 + -4 4 ,2,1 + 34,2,1* 

< 74 , 2,0 = -16/7X^2 ,_i32,2,2 - 16/7X^ 2,i32,2,-2 ■+• 8/7X2,2,-231,2,1 + 128 / 7 X 2,2 , -2-42,2,032,2,2 + 
16/7X2,2,-2-42,2,132,2,1 ~ 256/7X2,2,-2-42,2,232,2,0 ~ 64/7X2,2, -232, 2, o32, 2, 2 + 
48/7X2,2,-233,3,2 - 24/7X2,2,-lX2,2,o32,2,l + 48/7X2,2,-lX2,2,l32,2,0 + 
16/7X2,2,-lX 2 ,2,232,2,-l “ 8/7X2,2,-l3 2 ,2,-l32,2,2 + 12/7X2,2,-l32,2,o32,2,l + 

6X2,2,- l3s, 1,1 - 48/35X2, 2, -l3s, 3,1 - 24/7X2, 2, 0-42, 2,1 32, 2,-1 + 128/7X2,2,0-42,2,232,2,-2 + 
128/7X2,2,032,2,-232,2,2 ~ 24/7X2,2, 032, 2,^ -l32, 2,1 ~ 8/7X2,2,l32,2,-232,2,l + 

12/7X2, 2, l32,2,-l32, 2,0 “ 6X2,2,133,1,-1 + 48/35X2,2,1 3s, s,_i + 8/7 X2,: 2,232,2,-1 - 
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64/7X2,2,202,2,-202,2,0 - 48/7A3.3, 2-03,5, -2 + -^4,2,0 + -04,2,0. 

Ci , 2,-1 = - 88 / 21A ] 2 ,_ 1 52 ) 2,1 - 16/7X^2,0-02, 2,-1 + 256/21X2,2,-2X2,3-152,2,2 + 
64/21X2,2,-2-^2,2,0-02,2,1 ~ 32/21X2,2,-2-^2,2,1-02,2,0 “ 704/21X2, 2, -2^2, 2, 2^2, 2,-1 “ 

32/3X2, 2, -2^2, 2, -1-02, 2, 2 + 16/21X2,2,-2-02,2,0-02,2,1 + 8 X 2 , 2 , -2-03,1,1 +96/35X2,2,-203,3,1 + 
8/7X2,2,-102,2,0 + 16/7X2,2,— 1 X 2 , 2 , 002, 2,0 + 88/21X2,2,-1X2,2,102,2,-1 + 
64/3X2,2,-1X2,2,202,2,-2 + 352/21X2,2,- 102,2,-202,2,2 ~ 44/21X2,2,-102,2,-102,2,1 + 
2^2,2,-103,1,0 “ 48/35X2,2,-103,3,0 “ 32/21X2,2,0^2,2,102,2,-2 + 16/21X2,2,002,2,-202,2,1 “ 
8/7X2,2,002,2,-102,2,0 “ 4X2,2,003,1,-1 + 32/35X2,2,008,3,-1 + 44/21X2,2, 102,2,-1 ~ 
32/21X2,2,102,2,-202,2,0 ~ 128/21X2,2,202,2,-202,2,-1 “ 96/7X2,2,203,3,-3 + -44, 2,-1 + 04, 2 ,- 1 , 

£4, 2, -2 = 256 / 21 X 2 2, -202, 2, 2 — 8/21X3 3 _ 1 02,2,o ~ 64 / 21 X 3 , 2 , O 02 , 2, -2 + 32 / 21 X 2 , 2 ,- 202 , 2,0 — 
64/21X2,2,-2-42,2,-102,2,1 + 64/21X2,2,-2-42,2,002,2,0 ~ 16/3X2,2,-242,2,102,2,-1 ~ 
256/21X2,2,-2-42,2,202,2,-2 + 128/21X2,2,-202,2,-202,2,2 “ 88/21X2,2,-202,2,-102,2,1 + 
4X2,2,-203,1,0 + 24/35X2,2,-203,3,0 + 8/21X2, 2, -1X2,2, 002, 2,-1 + 176/21X2, 2, -1X2,2, 102, 2, -2 + 
8/3X2,2,-102,2,-202,2,1 “ 4/21X2,2,-102,2,-102,2,0 ~ 2X2,2,-103, 1,-1 ~ 24/35X2,2,-103,3,-1 + 
4/21X2,2,002,2,-1 “ 32/21X2,2,002,2,-202,2,0 + 8/7X2,2,003,3,-2 + 32/21X2,2,102,2,-202,2,-1 ~ 
24/7X2,2,103,3,-3 - 128/21X2,2,202,2,-2 + 44,2, -2 + 04,2, -2! 

<?4,0,0 = 44,0,0 + 04,0,0- 

This jato operation # is a noncommutative product, here expressed in the basis of the 
three-dimensional harmonic oscillator states. It has several properties that link the to physical 
properties of the optical elements with mathematical statements on selection rules [10]. But 
further, if we count the number of terms in the preceding gato operation in the symplectic bases, 
and a corresponding count in the monomial bases, we find some economy in the symplectic basis 


BASIS 


ORDER 3 ORDER 5 ORDER 7 


MONOMIAL 


SYMPLECTIC 12 


10 Economy in Aberration Calculations 

Two general linear (M) and aberration (A) group elements multiply through 

S{ A,M} £{B,N} = Q { A # D(MjJMN). (23) 


Now, the symplectic basis k Xm is block-diagonal under paraxial — linear — transfor- 
mations, with a matrix composed of the (analytically continued) Wigner D-mat rices: 
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The number of non-zero matrix elements in the two bases is 


BASIS 

ORDER 1 

ORDER 3 

ORDER 5 

ORDER 7 

k — ► oo 

MONOMIAL 

9 

36 

100 

225 

~ k*/ 4 

SYMPLECTIC 

9 

26 

58 

107 

~ 2k^ /Z 


The total number of operations necessary to 7 th aberration order, including matrix 
multiplication by the linear part is, for the two bases, 

monomial : 7680, symplectic : 3882. 

For this reason we conclude that the most efficient basis to carry through aberration compu- 
tations in axis-symmetric systems is the angular momentum basis of the harmonic oscillator 
states. 


11 Outlook 

This has been a quick revision of the state of the art in seventh order aberration calculations in 
axis-symmetric geometric-optics systems. We have seen that the state schemes of the harmonic 
oscillator provide order and symmetry in the classification of aberrations. In these conference 
proceedings we cannot go much further, so let me state that wave optics is the true objective of 
this quest. We can design and specify systems in geometric optics; once this is done we would 
like to predict the imaging behavior of such a system when light of a definite color is used. 
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